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Abstract . We generalize the previously proposed running vacuum energy model by 
including a term proportional to H, in addition to the existing term. We show that 
the added degree of freedom is very constrained if both low redshift and high redshift 
data are taken into account. Best- fit models are undistinguishable from ACDM at the 
present time, but could be distinguished in the future with very accurate data at both 
low and high redshifts. We stress the formal analogy at the phenomenological level of 
the running vacuum models with recently proposed dark energy models based on the 
holographic or entropic point of view, where a combination of H and term is also 
present. However those particular entropic formulations which do not have a constant 
term in the Priedmann equations are not viable. The presence of this term is necessary 
in order to allow for a transition from a decelerated to an accelerated expansion. In 
contrast, the running vacuum models, both the original and the generalized one in- 
troduced here contain this constant term in a more natural way. Finally, important 
conceptual issues common to all these models are emphasized. 



PACS numbers: 95.36.+X, 04.62.-Fv, ll.lO.Hi 
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1 Introduction 



The longstanding dark energy (DE) problem was originally presented in the form of the cosmo- 
logical constant (CC) problem [1]. Whichever way it is formulated, the CC problem appears as a 
tough polyhedric conundrum which involves many faces: not only the problem of understanding 
the tiny current value of the DE density pde in the context of quantum field theory (QFT) or 
string theory, but also the cosmic coincidence problem, i.e. why the density of matter pm is now 
so close to pde- Dynamical DE models are helpful in order to improve the situation. They can 
appear in different formulations of fundamental physics. Popular possibilities are, among others, 
quintessence and phantom energy in its various forms [2J , and scalar-tensor models [3| . Further- 
more, modified gravity is another very interesting option, which has been intensively explored in 
the recent literature, see e.g. [ll[5l[6]. 

But a class of cosmic accelerating models which we wish to explore in this paper is that of 
dynamical vacuum energy models. They have been proposed since long ago - see e.g. [3 El E] and 
references therein. Some of these "running" vacuum models are a possible clue for tackling one or 
more aspects of the CC problem. Despite the various phenomenological existing studies of time 
evolving vacuum models |10| . some of them are expected on more fundamental grounds, e.g. within 
the context of QFT in curved space-time [51 E]. In fact, it is difficult to conceive an expanding 
universe with a strictly constant value of the vacuum energy density p\ = A/{8ttG), namely one 
that has remained unchanged since the origin of time. It is much more natural to expect that 
the vacuum energy is a dynamical quantity as the universe itself, and thereby sensitive to time 
evolving functions such as the the Hubble rate H = H{t) or the scale factor a = a(t) = (1 + z)~^ 
(oo = !)• In these models, the need for scalars is obviated and nevertheless a phenomenologically 
viable description for the dynamical nature of the vacuum energy is achieved. Not only so, some of 
these models have been successfully tested against the latest cosmological data, see e.g. the recent 
studies im [T^ . Remarkably, some particular formulations of them have been used to improve both 
the cosmic coincidence problem [13] and the tough "old CC problem" , i.e. the fine tuning problem 
- see e.g. the recent attempts within the context of modified gravity [6]. 

More recently, Verlinde fl4J proposed that the gravitational field equations can be derived from 
the second law of thermodynamics in a way that would render the gravity force quite literally as a 
kind of "entropic force" (which is certainly not the case in e.g. Jacobson's |15) and Padmanabhan's 
approaches [16] . in which the entropic formulation is much more general). When Verlinde's entropic 
version is applied to cosmology, the DE does not exist anymore as such, but is mimicked in 
an effective way by the acceleration associated to the entropic force acting outwards the cosmic 
horizon. It is this particular formulation that can be called "entropic-force cosmology" which was 
first explored in Ref. [171 [18] and later on by various authors - see e.g. [191 |20l [21] . We emphasize 
that these models seem to lead to the same effective Friedmann equations as the aforementioned 
running vacuum models [71 [51 [S] with the notable difference that some of these entropic-force 
models, but not necessarily all of its versions, do not yield a constant term in their Friedmann 
equations. 

The plan of the paper is as follows. In Section 2 we review the running vacuum model followed 
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by a comparative discussion with the entropic- force models, and we emphasize the analogy at the 
level of the equations of motion. In Section 3 we present the background cosmology for these 
models. We show that the entropic-force cosmology appears as a particular case of the generalized 
running vacuum model. After comparing and fitting them to the data in Section [H we provide our 
discussion and final conclusions in Section [SI 



2 Running vacuum energy and entropic-force models 

As mentioned in the introduction, dynamical dark energy is an attractive possibility in order 
to explain certain aspects of the cosmological constant problem. In this section we review the 
idea of running vacuum energy, which was suggested in the literature long ago [22] , and we take 
opportunity to compare it with the more recent notion of entropic dark energy, specially some 
recent formulations of it [n\ I18j. They are formally similar but present also important differences 
which lead to significant phenomenological implications. The latter will be analyzed in subsequent 
sections. 

2.1 Running vacuum energy as dynamical dark energy 

The running vacuum energy in QFT in curved space-time derives from the renormalization group 
(RG) equation suggested in the literature for p\ - see [7] and references therein: 



d In ^2 (47r) 



n2/x' + 0(/), (2.1.1) 



where Mi are the masses of the particles contributing in the loops, and Bi,Ci, .. are dimensionless 
parameters. The equation (j2.1.ip gives the rate of change of the quantum effects on the CC as a 
function of the scale /u. Only the "soft-decoupling" terms of the form ~ p? remain in practice, 
as the Mf ones would trigger a too fast running of the cosmological term0. The approximate 
integrated form of ()2.1.ip is very simple: 

PA{H)=n^ + n2H^ + 0{H^), (2.1.2) 

where, following the aforesaid works, we have set p = H as the characteristic mass scale for 
FLRW-like universes, and we will neglect the (much) smaller higher order powers of H. Indeed, 
notice that only even powers are allowed by the general covariance, and hence no other if^^-terms 
beyond (not even H'^) can contribute significantly on the r.h.s. of equation ()2.1.2p at any stage 
of the cosmological history below the GUT scale Mx ^ Mp, so that we omit them. The additive 
constant term no in ()2.1.2p appears in a natural way in this framework upon integrating the RG 
equation. It will play a fundamental role in our discussion. Both no and ni become related by the 
boundary condition p\{Hq) = p\, which is to be satisfied by (j2.1.2p at present, Hq = H{to). As a 



^The main contribution to the running of pA clearly comes from the heaviest fields in a typical GUT near the 
Planck scale, i.e. those with masses Mi ~ Mx ^ Mp. See e.g. [8] for a specific scenario within this class of models, 
where the one-loop contribution for the Bi coefficients is explicitly given. 
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result these coefficients can be conveniently rewritten as follows: 

r^o = PA-^Ml, Hi , n2 = ^Ml, (2.1.3) 
where from (j2.1.ip we have defined the important dimensionless parameter 
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iy = — yB,^. (2.1.4) 

This parameter provides the main coefficient of the /3-function for the running of the vacuum 
energy. The coefficients Bi in (j2.1.4p can be computed from the quantum loop contributions of 
fields with masses Mj, and hence v is naturally expected to be non- vanishing and small {\v\ <C 1). 
For instance, for GUT fields with masses Mj near Mx ~ 10^^ GeV, a natural estimate lies in the 
approximate range v = 10~^ — 10~^ [8]. As a result we also expect a mild running of /OAi hence a 
dynamical DE framework which is healthfully close to the well tested concordance ACDM model 
in which pA is strictly constant. This particular situation is retrieved only for = 0, for which 
PA = Pa ^ times. However, there is no obvious reason for u to be strictly vanishing in QFT 
in curved space-time. Therefore in the general case we should have a time evolution law for the 
vacuum energy (j2.1.2p . whose leading contribution can be presented as follows: 

Pa(^) = pU ^ M| {H" - Hi) . (2.1.5) 

Substituting (j2.1.5p in the general acceleration law for a FLRW-like universe in the presence of a 
vacuum energy density pA; we find 

a AtkG 47r G 2 

- = — [Pm + 3pm - 2pA) = — (1 + 3uJm)Pm + Cq + U H , (2.1.6) 

a 6 6 

with 

^^0 = -i^Hl = I -i^Hl = Hi {nl - i.) . (2.1.7) 

Here ujm = Pm/ Pm is the equation of state (EoS) for a generic component of matter {ujm = and 
1/3 for non-relativistic and relativistic matter, respectively), and 0^ = A/ (SHq) is the cosmological 
CC parameter whose observational value is ~ 0.73. We note the presence of the constant term 
Cq (X uq. As warned before, this term will play an important role in our study. 

2.2 Entropic-force models and effective dark energy 

It is interesting that an effective dynamical dark energy component similar to the one derived in the 
previous section can also be motivated within the context of the entropic models. In a particular 
version of this framework, called the entropic-force models [13], the holographic screen is thought 
to induce a force F = T VS on a test particle near the screen, where T is the temperature of the 
screen and VS* is the change of entropy associated with the information contained in it (which 
involves a large number of d.o.f.). The screen is supposed to increase its entropy when the test 
particle approaches it. Therefore, VS* and the normal n on the screen (pointing towards the 
particle, located in the inner volume bounded by the screen) have opposite signs. Since the force 
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is directed towards the screen we have F = —TdS/dr, with dr the distance of the nearby particle 
to the screen. When appHed to cosmology |17]. the entropy of the Hubble horizon Rh = c/H 



the horizon increases by dr is simply dSn = 27r {Rh ks/lp) dr. Inserting it in the formula for the 
pressure exerted by the entropic force on the cosmological expansion, P = F/A = —{T/A) dSn/dr, 
and estimating that the horizon temperature is T = {h/kB){H /2t:) (proportional to the de Sitter 
temperature) pT], one finally obtains P = — (2/3) /JcC^, where pc = 3-ff^/(87rG) is the critical 
density. The minus sign in the pressure is of course the characteristic feature of the accelerated 
expansion in this entropic version. Apart from some coefficients that depend on the estimations 
made, and which are not essential for the argument, the basic result is that P oc —p'^^, where 
p'^^ ~ H'^ Mp (with Mp = G~^/^) is the quantity that plays the role of effective DE in this 
entropic model. This framework suggests that the entropic force leads to an effective DE density 
which is dynamical: it specifically evolves as the square of the Hubble rate[f|. By Priedmann's 
equation, it immediately follows that at the present time the value of p^^ would be predicted 
in the ballpark of the measured vacuum energy density: Pjy^ito) ~ H^Mp ~ p^. ~ 10"^"^ GeV^, 
where p\ = A/(87rG). 

Since the previous (entropic inspired) result is essentially a surface effect from the horizon, 
one may think of fully generalizing it by considering the gravitational action for space-times with 
boundaries [23]. This is achieved by adding the boundary action term Ip to the standard Einstein- 
Hilbert action, I eh, namely: 



Here h is the determinant of the metric hab on the boundary dM. , induced by the bulk metric g^y 
of A4, and y"" are the coordinates on dAi. Furthermore, K is the trace of the second fundamental 
form (or extrinsic curvature); if is the normal on the boundary, it can be written as K = V^n^. 
The complete action is / = I eh + Ib + Im, where Im represents the ordinary matter contribution. 
As a mere technicality, let us point out that the precise definition of the boundary term Ip should 
actually include an overall sign, which is plus or minus depending on whether the hypersurface 
dM. is space-like {n^n^ = +1) or time-like {n^n^ = —1), respectively. We exclude null surfaces 
for this consideration. Notice that the precise coefficient in front of the boundary integral Ip is 
chosen in such a way that the surface terms generated from the metric variation of Ieh are exactly 
canceled by the metric variation of Ip, provided the variation 5g^^ is performed in such a way 
that it vanishes on dAd, i.e. provided the induced metric hab on the boundary is held fixed. It 
follows that, in the presence of Ip, the standard form of Einstein's equations is preserved even if 
the space-time has boundaries. 

The authors of Ref. [17| presumably used the above interpretation of Ip as a way to generalize 
the entropic force argument given before, in the following way. As the surface terms emerging from 

^ Note that for the sake of better clarity, we keep h and c in this section, but natural units h = 1 = c for the rest. 
Apart from the running vacuum model of sect. 12.11 other frameworks involving the dynamical term ~ 
(treated as the full DE density or as a component of it) were suggested in ,24. and more recently in [25] , all of them 
involving the idea of vacuum fluctuations. 





(2.2.1) 
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the variation of I eh are canceled by 51 they assumed that if the total action would not contain 
Ib the contribution of the aforementioned surface terms to the field equations would be of the 
order of the effect induced on them by Ib, estimated as R times the prefactor l/(87rG) in Ib, i-e. 
(12ff^+6ij)/(87rG) - evaluated in the FLRW metric, in which H = a/a and H = dH/dt. However, 
since this is probably just a rough estimate of the effect, they finally proposed to generalize the 
corresponding acceleration equation for the scale factor in the form: 



However not all of the models considered in [17] are of this type Q The new ingredients are 
Ch and C^, which are certain (presumably small) dimensionless coefficients to be fitted to the 
observational data. Let us also mention that there can be higher order quantum corrections on the 
r.h.s. of Eq. (j2.2.2p - cf. [18]. We have neglected these effects for the present discussion because 
they have no impact for virtually any time in the history of the universe after inflation. This is 
in line with our approximation of ignoring the 0{H^) quantum corrections also in the running 
vacuum model discussed before. 

Let us point out that the field equations ()2.2.2p are not necessary derived from a fundamental 
action. Let us recall that in the most general entropic-holographic formulations, gravity is con- 
ceived as an emergent phenomenon [16] , and in this sense the gravitational field equations need not 
necessarily be deducible from a fundamental action at the present macroscopic level of description, 
even though the field equations themselves may provide a fully satisfactory account of all the basic 
phenomena known to date. From this point of view, the ultimate origin of gravity may lie in 
some fundamental degrees of freedom quite different from the metric variables, namely degrees 
of freedom which are completely unknown to us at present [TB]. If so, the field equations under 
discussion in this paper could just be effective field equations falling in this category and therefore 
no fundamental action to derive them would be needed. A detailed discussion on this point goes 
beyond the scope of the present work. 

The quantities H and appearing on the r.h.s. of (|2.2.2p are related through H = — (g+1) 
where q is the deceleration parameter. During some stages of the cosmic evolution when q is 
roughly constant, H and are approximately proportional. For example, (7 ~ 1 for the radiation 
dominated epoch, and g ~ 1/2 for the matter dominated epoch. Hence H ~ —2H^ deep in 
the radiation dominated era and H ~ —(3/2)//^ deep in the matter dominated epoch. When 
we compare the entropic formula (|2.2.2p with the corresponding equation (|2.1.6p in the running 
vacuum model, we see that deep in the matter dominated epoch we can set the correspondence 
I' Ch — 3Cjj/2 between the two models. However this is not valid at low redshifts when the 
universe goes over from matter domination to accelerated expansion. In this interval where SNIa 
data are located, q experiences a sharp model-dependent variation. Therefore the addition of the 
term C^H is a genuine extension of the original running vacuum energy model. 

In view of the close analogy between these models, in the next section we consider a gen- 
eralization of the running vacuum model with the inclusion of a term H together with the 




(2.2.2) 



one. 



*G. Smoot, private communication 
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3 Background solution of the cosmological field equations 



In this section we consider the solution of the cosmological field equations for both the generalized 
running vacuum model and the entropic-force model. We discuss in detail the underlying local 
conservation laws of matter and radiation in interaction with a dynamical vacuum energy compo- 
nent and we show that this leads to important conceptual issues. Finally, we emphasize the crucial 
importance of a constant term which rules out some of the entropic-force models lacking this term. 

3.1 The generic cosmological framework 

The cosmological equations of both the running vacuum models and the entropic-force models can 
be solved in a common framework. We will consider spatially flat Friedmann-Lemaitre- Robertson- 
Walker (FRLW) cosmologies 

ds"^ = dt^ - a^{t) dx^ . (3.1.1) 

Hence the (expansion) dynamics is fully encoded in the time evolution of the scale factor a{t). 
Instead of obeying the usual Priedmann equations of General Relativity, our models obey modified 
Friedmann equations, viz. 

= ^ + Co + ChH^ + C^H (3.1.2) 

i 

- = -^y,{p^ + ^Pi) + Co + CHH^ + CHH , (3.1.3) 
a 3 ^—^ 

where the remaining sum is over the matter components only. For realistic cosmologies, we take 
as usual two components, namely nonrelativistic (dust-like) matter with pm = Wm Pm = and 
radiation with = Wr Pr = \ Pr- These equations can be viewed formally as resulting from the 
presence of a time-dependent component pA(t) = pA{H{t), H{t)) satisfying 

Pa{H, H) = ^ (Co + ChH^ + Cjjii) = -pa{H, H) . (3.1.4) 

We will call this dynamical component a "generalized running vacuum energy" (GRVE) density 
since its EoS satisfies w\ = pa/ PA = —1 as in the case of a strictly constant vacuum energy[f| . 
By the same token we will call the class of these models with Cq ^ the "generalized running 
vacuum models". In the particular case = we recover the original running vacuum model 
discussed in sect. 12.11 Formally, the generalization of the model being proposed here implies that 
the scale /i^ in Eq. ()2.1.ip , which is to be eventually associated with a physical quantity according 
to the RG procedure, should in general be a linear combination of and H rather than just 
the component, as these two terms represent independent d.o.f. with the same dimension. 
Finally, let us emphasize that the particular case Co = is not to be included within the class 
of GRVE models because the integration of the RG equation (j2.1.ip always involves an additive 

^Let us notice that the recent work [26], which extends the discussion of a model first suggested in Appendix C 
of [11], contains a linear term in H rather than our H term. As mentioned in sect. 12. 11 odd powers of H cannot 
emerge from a covariant effective action, and in this sense these models are more phenomenological than the class 
of GRVE models (|3.1.4I) presented here. 
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term leading to Cq ^ 0. The case Cq = seems to appear in some of the entropic-force models [T7] 
briefly addressed in the previous subsection. While this setting can be derived as a particular case 
of our general analysis of the system (j3.1.2p - (|3.1.3p . we stress that Cq = leads to a qualitatively 
new situation, which we will comment in subsequent sections and that is not expected from the 
conceptual point of view of the running vacuum model framework. With these provisos in mind 
we are now going to solve the background cosmology of the entire class of models (I3.1.2p - (|3.1.3p . 

3.2 Discussion of the local conservation laws 

Once the metric (j3.1.ip is given, a comoving perfect fluid with energy-momentum tensor T^u = 
{Pi + Pi) UfiUu — Pi Ofiu will satisfy the conservation equation pi = —3H{pi + Pi), with pi and pt 
appearing in ()3.1.2p - ()3.1.3p provided pA is constant. This applies of course for non-flat FLRW 
universes as well. However, if the vacuum energy density is a time-dependent component it cannot 
have this energy-momentum tensor. Actually, as we will see now the same applies for the other 
components appearing in equations (|3.1.2p - (|3.1.3p . 

However, from the system (|3.1.2p - (|3.1.3p we get the coupled conservation equation 

Pm + Pr+ PA = -^Hpm-4.Hpr. (3.2.1) 

This equation is indeed a first integral of that system. Clearly, for PAj^ none of the components 
can satisfy the standard conservation equation as emphasized above|j. 

It is easy to see that this does not depend on the particular choice (|3.1.4p . Let us assume that 
the variable vacuum has an energy density ^ pA{t) = Cr H"^ and that it is a perfect fluid with 
an a priori undeflned equation of state. Then it follows from its (assumed) conservation equation 
that its equation of state parameter wa satisfies —jjr = |(1 + ^^a)- So if we consider a simple 
(flat) universe containing also dust, the only consistent way is to have = Ch = 1- So we end 
with no dust at all and the flrst Priedmann equation reduces to an equality (while wa remains 
of course undeflned)! We stress that these conceptual issues are also true for other models with 
similar effective Priedmann equations, for example models inspired by the holographic principle or 
models based on the entropic-force principle. Ultimately these properties arise from the absence 
(in general) of a formulation at the level of the action which is still an open issue. 

As mentioned in the previous section, in a pragmatic approach we assume the validity of 
(j3.1.2p - (j3.1.3p without explicitly deducing them from an underlying action - see, however, |8j for 
a speciflc framework along these lines. After some calculations the coupled conservation equation 
(|3XT]) reads 

Pm + Pr - -^Cjj (pm + -Prj = -3H{1 - Ch) ipm. + -^Prj ■ (3.2.2) 



^The non-conservation of matter in the presence of running vacuum energy has recently been proposed as a 
possible link between the dynamical DE and the increasing evidence for a possible variation of the fundamental 
constants and scales in Nature, as e.g. the QCD scale - see [27] • However, a running vacuum energy of the form 
(|2.1.5|l can be made compatible with matter conservation if one allows G to slowly evolve with time, see [8! for a 
concrete scenario connected with an action functional. For the present GRVE framework, though, G is assumed to 
be strictly constant. 



8 



We note that the model does not yield a conservation equation for each component separately, 
seemingly overlooked in[T7]. For this we would need to specify the action of the model and to 
find the corresponding energy-momentum tensors. Moreover, the system is not fully defined by 
eqs. (|3.1.2p - (j3.1.3p . Indeed, any solution of the following set of equations 

1-Ch 

Pm = -3H- 3— Pm + Q (3.2.3) 

Pr = -4g ^_^J^ Pr-Q , (3.2.4) 

for arbitrary function Q = Q{t), will be a solution of equation (j3.2.2p . In the matter dominated 
(pr ~ 0) and radiation dominated (pm ~ 0) stages we must have from p.2.2p that Q ^ 0. The 
simplest version of this model is to assume that eq. ()3.2.2p reduces at all times to a set of decoupled 
equations with Q = at all times. We conjecture that this is perhaps the only way to introduce 
consistently an arbitrary number of species. 

It will be convenient to introduce the following notations 



V = Ch (3.2.5) 

3 
-( 

2 



a = -Cjj , (3.2.6) 



as well as the important quantities 



U = r. (3.2.7) 

\ — a 

1 — I' , ^ 

& = — ■ (3-2-8) 

The motivation for the relabeling (j3.2.5p is simply because for = the GRVE model boils 
down to the original running vacuum model discussed in sect. 12. 11 and then Ch exactly reduces 
to the parameter u defined in that section. Using these definitions, equations (I3.2.3P and (|3.2.4p 
(setting Q = 0) can be written as follows 

Pm = -3H Pm (3.2.9) 
p^ = -AH^r pr , (3.2.10) 

for which it is straightforward to obtain the corresponding solutions (setting uq = 1): 

Pm = a-^?'" = (1 + z)3€'" (3.2.11) 
p^ = pO a-^«'- =p° (l + z)^^^ . (3.2.12) 

Note that these decoupled solutions reduce automatically to the behavior of dustlike matter during 
matter domination (p.,. ~ 0) and to the radiation component during radiation domination (pm ~ 0). 
They take the standard form for (,m = ^ and = 1- 
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3.3 Determining the time evolving vacuum energy and the Hubble function 

The equations (j3.2.9p and (j3.2.10p are decoupled, there is no transfer of energy between the two 
components. However there is a transfer of energy between the running vacuum energy p\ and 
these components. We find for the evolution of p\ 

PA = 3H{U -1) Pm + ^H{ir -I) Pr . (3.3.1) 

When Cm = = 1 the standard behavior of matter and radiation is recovered and then = 
A/{8ttG) reduces to a genuine cosmological constant with A = 3 Cq. The transfer of energy 
between the matter components and the GRVE is the physical reason for the particular scaling 
behaviors ()3.2.1ip and ()3.2.12p . which obviously depart from the standard expectations pm ci~^ 
and Pr in the ACDM owing to the non-vanishing values of the parameters u and a. We 

have here an effective interacting dark energy model. Consistency enforces to have an interaction 
between the running vacuum energy and all other components. 

The addition of the term H = H introduces an important change compared to the 
original running vacuum model discussed in sect. 12. II - in which = - but the new degree of 
freedom is severely constrained by observations. The reason is that the model cannot depart too 
much from the ACDM values Cm = = li with 

Cr = Cm ~ 4 • (3.3.2) 

1- gO 

The only way to satisfy both constraints Cm ^ 1 and Cr ~ 1 is to have 

|i/|<l, |a|<l ^ Cm~l, Cr~l (3.3.3) 

Of course we have = = ^ when both parameters and a vanish. Note also that the condition 
Cr ~ 1 is crucial for the viability of our model at high redshifts (e.g. when fitting the model against 
CMB data). If we use only constraints at very low redshifts, models with Cm ~ 1 (i-e. <SC 1) but 
not necessarily satisfying |a| ^ 1 will fare well in this domain, though the model would actually 
be unviable taking into account its behavior at high redshifts. 

Compared to the old running vacuum energy model (a = 0), the generalized model offers more 
possibilities to depart from standard cosmology: 

• One can have Cr = Cm ~ !> in which case the model just reduces to the original running 
model (a = or = 0, and u ^ 0). Both radiation and dust scale in a non-standard way 
but their departure from standard behavior is not independent and depends on one single 
parameter ly, specifically: pm ~ a^^^^^'^^ and pr ~ a^^^^^'^\ 

• Cr = 1 and Cm ^ 1, in which case radiation behaves in the standard way but dust does 
not. This occurs when = |a 7^ 0. This case exists only in the generalized model a 7^ 0. 
Departure from standard cosmology occurs already at low redshifts. 

• Cm = 1 and Cr ~ 1, now dust scales in standard way but radiation does not. This corresponds 
to v = a ^ 0. This case can mimic standard cosmology at low redshifts but it is strongly 
constrained when high redshift data are considered; 
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• '?m ~ li ~ 1 Eiiid Cr- Here the deviation of the non-relativistic component is different 
from the relativistic one, and hence this provides an extension of the first case discussed 
above which is only possible within the GRVE model. 

The three last cases above are only possible in the generalized model a ^ 0. However these 
additional possibilities are strongly constrained by observations on both low and high redshifts. 
We will see in particular in sect. U] that the strong constraint on a in the regime (j3.3.3p is similar 
to that one found for in the original running model - see the recent analyses [ll^ [12] . 

Equation (13.3. ip is easily recast in the form 

^ = HU-l)—+^{Cr-l)^ (3.3.4) 
da a a 

which is easily integrated using the solutions (j3.2.1ip and (|3.2.12p and the explicit form of the 
vacuum energy as a function of the scale factor can be expressed as follows: 

Pa(«) =pi + p'L (C' - 1) {a-'^- - l) + Pr (C' - 1) («-'^^ - l) • (3.3.5) 

The Hubble function can now be constructed from the matter components (|3.2.9p - (|3.2.10p and the 
vacuum energy (j3.3.5p : 



= [Prn + Pr + Pa] = ^ 



qO ^ I qO ^ , qO P± 



(3.3.6) 



Pm Pr Pa . 

Introducing the normalized Hubble rate in terms of the redshift, E{z) = H{z)/Hq, we find: 

E\z) = ^ (1 + zf^- + + z)^^^ + , (3.3.7) 

where we have used the standard definition = Pi/pc, with pc = 3H'^ /{SirG), satisfying the 
constraint 

nA + Qrn + ^r = 1 , (3.3.8) 

at all times. 

Note that the boundary condition E{z = 0) = 1 in (|3.3.7p leads to the equality 

Hq^Cq = nl-Au, (3.3.9) 

where we have defined Aiy = v — i?, with P = a + (4/3) a $7^^. The P parameter is characteristic 
of the extension of the original running vacuum model into the GRVE model and is closely related 
to u. Indeed u gauges the size of the new i7-effect in terms of at the present time since it 
satisfies the relation C^Hq = —uHq, which can be compared to Ch Hq = Hq in the original 
running model. To confirm that relation let us write the current value of the C^H term in the 
starting equations (j3.1.2p - (l3.1.3p as follows 

CjjH^ = -C^{q, + l)Hl = -{^-nl + 2^i^ C^Hl (3.3.10) 

Thus we find u = i?^ + v^, where Vm = (3/2) O^C^ = "f^m and Dr = 2QP^Cjj = (4/3)aJ70 
represent the non-relativistic and relativistic matter contributions respectively. 
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After having determined the explicit relation between v and the other parameters, we see from 
Eq. (j3.3.9p that Co becomes also explicitly determined as follows: 



Cn 



0° 



a 



(3.3.11) 



Notice that for a = (or = 0) it boils down to the corresponding expression (I2.1.7P for 
the original running vacuum model. On the other hand Eq. ()3.3.1ip tells us another interesting 
feature, to wit: models with Co = cannot have the two parameters v and a (equivalently Ch 
and Cjj) simultaneously small, i.e. it is impossible to satisfy the relations (|3.3.3p . unless = 
- which is of course unacceptable. In particular, entropic-force models [TTj cannot have Ch and 
C^ simultaneously small, otherwise they would contradict the measured value of the cosmological 
term: il^ ~ 0.73. Even if we would accept that at least one of the parameters Ch and C^^ is not 
small, the resulting model would be contrived as it would entail a non-trivial modification of the 
standard ACDM cosmology. Actually in the next section we will encounter a related difficulty, 
which is perhaps the biggest stumbling block to the Cq = models. 



3.4 Crucial distinction between some entropic-force and GRVE models 

In the previous subsections we have solved in detailed the full class of cosmological models based on 
the set of generalized FLRW equations (j3.1.2l) - (l3.1.3p . In particular, we have assumed arbitrary 
values for the parameters Co, Ch = v and C^ (x a. However we expect from observational 
constraints that the last two ones are sufficiently small ~ cf. Eq. (j3.3.3p - in order for the generalized 
models not to depart too much from the standard scaling laws of matter and radiation. 

We turn now our attention specifically to the additive parameter Co- If the other two param- 
eters {Ch and C^) have to be small, this is not the case for Cq and as we will see now it cannot 
vanish. While the running vacuum energy models have a nonvanishing Co, this is not the case for 
some entropic force models. 

Indeed, successful models must be able to produce an accelerated expansion at very low red- 
shifts. To start with let us analyze the situation Cq = 0. It is easy to derive from the expression 
for I that accelerated expansion is obtained both in the matter and radiation-dominated stages if 
the following condition is satisfied (with Cq = 0) 

z.-^a>^ ^ 2ir<l. (3.4.1) 

In the matter-dominated era a slightly weaker condition is required 

u-la>^ ^ ^em<l. (3.4.2) 

We see that these conditions are redshift independent. Therefore, if we have an accelerated expan- 
sion rate at very low redshifts, we will have it at least during all of the matter-dominated stage. 
This leads obviously to an unviable cosmology putting aside the fact that the corresponding scaling 
behaviors are completely unviable observationally. 
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We can recover these results solving for the time dependence of H(t) and a{t). The following 
equation holds during matter domination 



(3.4.3) 



When Co = the condition (13.4.2P for accelerated expansion is clearly recovered from (I3.4.3p . In 
the general case Co 7^ 0, equation (j3.4.3p can be solved to yield 



H{t) 
a{t) 



A coth 



D sinh^Jn 



where we have used (13.3.90 and we have set 



2 Cm A t 



A = Hn 



(3.4.4) 
(3.4.5) 



Qa,o ~ J \ \ — V 

Returning to the case with vanishing Co, the solution to (I3.4.3P reads 

2 



a{t) oc t-*?" 



3i^^ 



1 t- 



(3.4.6) 



(3.4.7) 



which shows again that (j3.4.2p leads to accelerated expansion. Hence we conclude that models 
with Cq = cannot describe an expanding universe undergoing a transition from decelerated to 
accelerated expansion. 



3.5 Observational interpretation 

For an "Einsteinian" interpretation of the generalized running vacuum models, i.e. those repre- 
sented by equations (|3.3.6p - (|3.3.7p with Cq / 0, let us write 



+ ~p,] + ^ = Hi 



qO Prn_ , qO Pr 

~o "^"i" p,0 
rm rr . 



A 

+ 3' 



with the obvious identifications 



Pm 

Pm — 7 , 



Pr A 



Cn 



Cm 



P'c 



Cr 



Q ^ ^m Pr 

^' ~ ^ ~ P'c 



(3.5.1) 



(3.5.2) 



Observationally there is no reason to distinguish between the matter or radiation energy density 
appearing in the starting equation (j3.1.2p and that part contained in p\. Hence it is natural to 
identify the observed value with and similarly ^2^.^^^^ with Cl^ . We still have the standard 

equality valid at all times 

nrn + ^r + ^A = l , (3.5.3) 

with = gA^. 

Even recast in the form p.5.ip we should remember that the energy densities pm and pr, obey 
the nonstandard scaling laws (j3.2.1ip . resp. (|3.2.12p . Interestingly, in this "Einsteinian" interpre- 
tation, our model reduces to a model with a genuine cosmological constant A and nonstandard 
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evolution of dust and radiation. In the generalized running vacuum energy model, the departure 
from standard behavior of dust and radiation are independent from each other. 
In this model, the redshift at equality z^q is given by 

(l + Ze,)'^'-3«'" = 2^ (3.5.4) 

The variation of Z(,q constrains the quantity 4^^ — 3^m) a constraint that will be satisfied by our 
best-fit models found in next section. In view of (j3.3.7p we expect further very tight constraints 
on ^j. itself deep in the radiation dominated era. 

In a first conservative approach we would like to keep a standard thermal history. Even if (l^ 
assumes the value for Q.^ required by standard Big Bang Nucleosynthesis (BBN), and assuming 
that cosmic temperature scales in the standard way, the expansion rate at the BBN epoch will get 
changed by a non-standard amount due to the scaling law (j3.2.12p . Inserting numbers this finally 
yields the conservative constraint j^^ — 1| < 10~^ because the expansion rate is severely constrained 
and cannot vary too much at the time of BBN (see e.g. |28j). 

As is a free parameter, in practice we wish to explore scenarios satisfying 

^,. = 1. (3.5.5) 

This choice means that the model parameters u and a are no longer independent, and from p.2.8p 
we see that we must have a = 2>v/A. This ensures that the standard thermal history is recovered. 
Indeed, with the choice (j3.5.5p . the radiation dominated stage in our models is essentially similar to 
the standard radiation dominated stage. We have in particular that the temperature of thermalized 
relativistic species scales consistently in the standard way. This is in particular true for the Cosmic 
Microwave Background (CMB) temperature. 

We have derived all equations necessary in order to constrain with observations the class of 
generalized FLRW models (j3.1.2p - (j3.1.3p satisfying (j3.5.5p . This we do in the next section. 

4 Fitting the models to the observational data 

In the following we present some details of the statistical method and on the observational samples 
and data statistical analysis that will be adopted to constrain the models presented in the previous 
sections. We shall extract our fit from the combined data on type la supernovae (SNIa), the data 
on the Baryonic Acoustic Oscillations (BAOs), and the shift parameter of the Cosmic Microwave 
Background (CMB). Note that in the case of the BAO analysis we have to modify it appropriately 
in order to incorporate some specific features of the present models. 
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4.1 The global fit to SNIa, BAOs and CMB 

First of all, we use the Union 2 set of 557 type la supernovae of Amanullah et al. [230. The 
corresponding x'^-function to be minimized is: 

557 

XSNIa(P) = X] 

4 = 1 

where Zi is the observed redshift for each data point. The fitted quantity [i is the distance modulus, 
defined as jj, = m — M = Slogd^ + 25, in which p) is the luminosity distance: 

dz' 

dUz,p) = {l + z) (4.1.2) 

Here p a vector containing the cosmological parameters of our model that we wish to fit for. In 
our case one possibility would be to take e.g. p = In equation M.l.ip . the theoretically 

calculated distance modulus ^th for each point follows from ()4.1.2p . in which the Hubble function 
H{z) = HqE{z) is given by ()3.3.7p for the generic model under consideration. Finally, /iobs(-Zj) 
and (Tj stand for the measured distance modulus and the corresponding la uncertainty for each 
SNIa data point, respectively. The previous formula (j4.1.2p for the luminosity distance applies only 
for spatially flat universes, which we are assuming throughout. Note that since only the relative 
distances of the SNIa are accurate and not their absolute local calibration, we always marginalize 
with respect to the internally derived Hubble constant (for methods that do not need to a priori 
marginalize over the internally estimated Hubble constant, see for example [30^ I31j). In the case 
of the Union2 SNIa data the internally derived Hubble constant is Hq ~ 70Km/s/Mpc which is in 
agreement to that of WMAP7 |38j Hq = 70.4Km/s/Mpc used in the present study. 

In addition to the SNIa data, we also consider the BAO scale produced in the last scattering 
surface by the competition between the pressure of the coupled baryon-photon fluid and gravity. 
The resulting acoustic waves leave (in the course of the evolution) an overdensity signature at 
certain length scales of the matter distribution. Evidence of this excess has been found in the 
clustering properties of the SDSS galaxies (see [32], [531131]) and it provides a "standard ruler" 
that we can employ to constrain dark energy models. In this work we use the results of Percival et 
al. [33], rs{zd) / DY{zi,) = 0.1390 ± 0.0037. Note that rs{zd) is the comoving sound horizon size at 
the baryon drag epoch [35] (i.e. the epoch at which baryons are released from the Compton drag 
of photons), and ~ 0(10^) is the corresponding redshift of that epoch, closely related to that 
of last scattering- the precise expression being given by the fltting formula of [35]. Finally, Dy{z) 
is the effective distance measure [32] and z^ = 0.275. Of course, the quantities {rsjDy) can be 
defined analytically. In particular, rs[zd) is given by the comoving distance that light can travel 
prior to redshift z^: 

rs{zd) 

where = (1 + -Zd)^^, and 

Cs{a) = 

^Note that the data can be found in: http: / / supernova.lbl.gov /Union / [ 



(4.1.1) 



a Jq a?H{a) 



(4.1.3) 



5p^ 



1/2 



6p^ + Spb J y/3lTTTZ(aj) 



(4.1.4) 



15 



is the sound speed in the baryon-photon plasma. Here we assume adiabatic perturbations and 
we have used 6ph = and dp.y = {l/3)6p-y, and defined TZ{a) = 6pb/5p-y. If the scahng laws for 
non-relativistic matter and radiation would be those of the standard model, we would have TZ{a) = 
3pb/4:py, which can be finally cast as a linear function of the scale factor: Tl{a) = (30.^/40,^^ a, 
where ~ 0.02263 and Cl^ h? ~ 2.47 x 10~^ are the current values of the normalized baryon 

and photon densities. However, our scaling laws for non-relativistic matter and radiation are given 
by equations ()3.2.1ip and ()3.2.12p . As a result, the sound speed velocity in the plasma gets a 
correction with respect to the standard result (oq = 1): 

^lUM^ 3 l-4a/3^ 4e.-3U . (4.1.5) 
^ ' A ir P^{a) 4 1 - a > 

Of course for v = Q and a = (^m = = 1) the previous equation becomes again a linear function 
of the scale factor, and it exactly reduces to the standard result. 

The remaining ingredients of the BAO analysis are as in the standard case, in particular the 
effective distance is (see [32]): 



Dviz) s 



1/3 

(4.1.6) 



where Da{z) = (1 + z) '^dL{z,p) is the angular diameter distance. Therefore, the corresponding 
Xbao function is simply written as: 

^^^^(p) -0.1390l^ 



xbao(p) = om37^ ■ ^ ' 

Furthermore, a very accurate and deep geometrical probe of dark energy is the angular scale 
of the sound horizon at the last scattering surface, as encoded in the location ij^ of the first peak 
of the Cosmic Microwave Background (CMB) temperature perturbation spectrum. This probe is 
described by the CMB shift parameter [361 EZ]; defined as: 

The measured shift parameter according to the WMAP 7-years data [38] is i? = 1.726 ± 0.018 at 
the redshift of the last scattering surface: zis = 1091.36. In this case, the x^-function is given by: 

2 , , [i?(p) - 1.726]2 ^ 
Xcmb(p) = ' Q Q^g2 • 

For a detailed discussion of the shift parameter as a cosmological probe, see e.g. [39]. Let us 
emphasize that when dealing with the CMB shift parameter we have to include both the matter 
and radiation terms in the total normalized matter density entering the E{z) function in (j4.1.8p . 
given explicitly by (|3.3.7p . Indeed, the radiation contribution reads VL^ = {\ + Q.221Ny) ilJJ, with 
Ny the number of neutrino species. Therefore, at zis = 1091.36, and including three light neutrino 
species, the radiation contribution amounts to ~ 24% of the total energy density associated to 
matter, which is not negligible. We use h = 0.704 in our analysis. 
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Our statistical analysis, due to its simplicity, has been used extensively in the literature in order 
to constrain the dark energy models (see for example [4:0\ \3T\ BT] and references therein) . We would 
like to point that a more general statistical presentation would require the covariances of BAO 
and CMB shift parameter. We have checked our statistical results using the latter covariances 
and our results remain the same as they should. Note that the corresponding covariances can be 
found in Percival et al. [33] and in Komatsu et al. |38j respectively. Finally, as emphasized before 
eq.(3.1.1) we restrict our analysis to spatially fiat spaces. This seems justified in view of the tight 
constraints on fifc g and is sufficient for our purposes. 

4.2 Numerical results 

Since we perform an overall fit of the SNIa+BAO+CMB data, it is important to take into account 
the contribution of both non-relativistic matter and radiation. 

• For the concordance ACDM cosmology, we simply have =const. and 

Pm{z) = p'^{l + zf, /5,(z) = pO (1 + ^)4. (4.2.1) 

• Concerning the generalized running vacuum energy model (|3.1.2p - (j3.1.2p we have explicitly 
given the corresponding density formulae in sect. El Let us recall that Cq = is a very 
particular case that we exclude from the class of the GRVE models. We have shown in the 
previous section that this case is not viable observationally and therefore we will not consider 
it any further for the phenomenological analysis. Therefore, from now on we assume that 
Co 7^ and focus on fitting the parameters of this model to the SNIa+BAO+CMB data. In 
particular, we already know that v and a have to be small - see Eq. ()3.3.3p - but only the 
direct confrontation of the model with the data will tell us about their possible maximum 
size. In practice, considering models satisfying (j3.5.5p . it will be convenient to define the 
effective parameter 

u^s = iy-a = -u, (4.2.2) 

and use u^g as fitting parameter, together with Q^. We can check this explicitly by expanding 
and linearly in <C 1 and |a| <C 1, together with some coefficients in (I3.3.5p . and 
using the definition (14.2.2P : 

^m^l-l^eS, = 1 . (4.2.3) 

As a result, for the energy densities we find: 

Pm{z) = pliz) (1 + z)3(i--=ff) , p,(z) = (1 + zf . (4.2.4) 

Similarly the corresponding normalized Hubble fiow squared reads 

E\z) = 0^ (1 + z)3(i"-cff) + 0° (1 + zf + nl . (4.2.5) 

Finally, let us mention that within the same approximation we can write the BAO ratio 
(|4.1.5p entering the modified sound speed of the baryon-photon plasma as follows: 

3 

7^(a) = - (1 - z^eff) ^ ai+='^- . (4.2.6) 
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For = and = (hence i^eflf = 0) it clearly reduces to the standard result mentioned in 
the previous subsection. 

These formulae confirm our contention that we can fully reexpress all the background for- 
mulae in terms of the effective fitting vector 

Peff= (j^^,l/eff) . (4.2.7) 

We also see from the previous formulae that i^eff is the single effective parameter that controls 
the deviations of the GRVE model with respect to the ACDM model in the low z region (when 
radiation can be neglected). It is only in the high redshift region where the model is sensitive 
to independent contributions from (equivalently, from v or a). Notice that this feature 
could be used, in principle, to distinguish between the two sorts of running models, i.e. the 
original one (which we reviewed briefly in sect. 12. ip and the generalized running vacuum 
model under discussion in this paper. At low z the two kinds of models are indistinguishable 
because they both depend on a single parameter, v and v^^ respectively. 

Let us next proceed with the numerical fit analysis. In order to place tighter constraints on 
the corresponding parameter space of our model, the probes described above must be combined 
through a joint likelihood analysifl given by the product of the individual likelihoods according 
to: 

Aot(p) = >CsNIa X L^KO X -CcMB , (4.2.8) 

Since likelihoods are defined as Lj oc exp ^— x^/2^ , it translates into an addition for the joint 
function: 

X?ot(p) = XsNIa + XbAO + XCMB • (4-2.9) 

In our minimization procedure, for the vacuum models (running and concordance ACDM) 
we use the following range and steps for the fitting parameters: G [0.01,1] in steps of 0.001 
and I'efr € [—0.02,0.02] in steps of 10~^. The numerical results that we obtain are the following. 
In the case of the generalized running vacuum model the overall likelihood function peaks at 
= 0.274 ± 0.011, z^efi = -0.00133 ± 0.0028 (or v = 4i^eS ^ -0.00532, a = Si^eS ^ -0.004) 
with Xtot(^m> '^cff ) — 542.93 for 556 degrees of freedonJ§. In Fig. [1] we present the la, 2a and 
3a confidence levels in the plane. In particular, the left panel in that figure shows the 

individual likelihood contours, with the SNIa-based results indicated by thick solid lines, the BAO 
results by dotted-red lines and those based on the CMB shift parameter by dashed-green lines. 
Using the SNIa data alone it is evident that although the parameter is tightly constrained 
(~ 0.27), the i^eff parameter remains completely unconstrained. As can be seen in the right plots 
of Figure 1, the above degeneracy is broken when using the joint likelihood analysis, involving all 

^Likelihoods are normalized to their maximum values. In the present analysis we always report la uncertainties 
on the fitted parameters. Note also that the total number of data points used here is Ntot = 559, while the associated 
degrees of freedom is: d.o.f = Ntot — n^t — 1, where riflt is the model-dependent number of fitted parameters. 

^Note that in [12] the original running vacuum model was used - see Eqs. (|2.1.5p and p.l.6p in which a is 
strictly equal to zero - and the Constitution set of 397 SNIa data |42) . We would like to mention here that those 
results for v are in agreement with the current results for foff within la uncertainties. 
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Figure 1: Likelihood contours (for — 2In£/£max equal to 2.30, 6.16 and 11.81, corresponding to la, 2a and 3(t 
confidence levels) in the (fjJij, j/efi) plane for the generalized running vacuum model (lITSll (Co / or ni 7^ 0). 
For the CMB analysis we include also the radiation component as indicated in (|4.2.4|l - (|4.2.5|) . The left panel shows 
the contours based on the SNIa data (thick solid black lines), BAOs (dotted- red lines) and CMB shift parameter 
(dashed-green lines). In the right panel we show the corresponding contours based on the joint statistical analysis 
of the SNIa+BAO+CMB data. 

the cosmological data. Finally, in the case of the concordance ACDM cosmology (i^eff = a = 0) we 
find = 0.274 ± 0.01 with xloti^m) / d-o-f - 543.18/558. 

Overall we see that the departure of the GRVE model with respect to the ACDM is extremely 
small and cannot be detected at present. 

5 Discussion and conclusions 

In this paper we have generalized the running vacuum energy models and we have solved the 
corresponding background cosmology. The generalized running vacuum model ()3.1.2p - ()3.1.4p with 
Co 7^ is able to pass the SNIa+BAO+CMB data constraints with a statistical significance 
comparable to that of the concordance ACDM model which is a limiting case of the model {^r = 

= 1, or equivalently a = v = {}). Although the best-fit models are currently indistinguishable 
from ACDM we expect that future very accurate data on both low and high redshifts could help 
to distinguish these models from the standard cosmology. 

Some conceptual issues pertaining to these models were also addressed which are related to 
the peculiar conservation laws derived in Section 3. We stress that these issues hold as well 
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for other models with analogous effective Friedmann equations, models inspired either by the 
holographic or the entropic-force principle. We have further emphasized that the presence of a 
non vanishing additive constant Co is crucial since otherwise the cosmology does not allow for a 
transition between decelerated and accelerated expansion. That was actually noticed in previous 
entropic-force studies [2D]0. In contrast, the class of the running vacuum models, both the 
generalized one (GRVE) presented here and the original one (which existed in the literature since 
long ago - see [7] and references therein) do not suffer from this problem because Cq is naturally 
expected to be non-vanishing as a result of integrating the corresponding RG equation. Therefore, 
despite the formal analogies between these two sorts of models, the running vacuum models are 
naturally well positioned for a correct phenomenological description of our cosmos. 

From the point of view of the running vacuum models, the current Universe appears as FLRW- 
like with a genuine cosmological constant while dust and radiation evolve in a nonstandard way, 
in the sense that they follow scaling laws that deviate slightly from their behavior in ACDM. In 
contrast to the old running vacuum energy model, the generalized one introduced in this paper 
allows for an independent departure from the standard behavior of both components. We have 
used this freedom and we have explored models satisfying = 1 thereby ensuring that relativistic 
matter obeys the standard behavior. In this way potential difficulties related to the radiation 
dominated era are essentially avoided. While the other parameter ~ 1 — fgfj remains free, it 
can be efficiently constrained using CMB data. It is constrained by observations at a similar level 
as the single parameter of the original running model, i.e. they are both presently allowed up to 
0(10^^) at most (in absolute value). 

This order of magnitude size is consistent with the theoretical expectations on these coefficients, 
interpreted as one- loop /3- functions of the running cosmological constant. The mild variation 
induced on the CC term by these coefficients is responsible for the dynamical character of the 
vacuum energy, which is of course the reason why these models have a chance to improve the 
situation with the ACDM without giving up its phenomenological success. Such time variation is 
foreseen on general QFT grounds and it provides a possible formulation of an effective dynamical 
dark energy, which in some cases can help curing the cosmic coincidence problem [33] and other 
related problems. 

To summarize, the running vacuum models offer a challenging phenomenologically consistent 
description of a universe with presently accelerated expansion. The dynamical A could be under- 
stood in the context of QFT in curved space-time. Such potential connection with fundamental 
physics could help to conceive the origin of a dynamical A term in QFT and eventually provide an 
explanation for the tough cosmological constant problem. 

Acknowledgment s 

JS has been supported in part by MEG and FEDER under PA2010-20807, by the Spanish 
program GPAN GSD2007-00042 and by 2009SGR502 Generalitat de Gatalunya. SB thanks the 
Dept. EGM of the Univ. de Barcelona for the hospitality, and the financial support from the 
Spanish Ministerio de Education, within the project SAB2010-0118. 

^"For other problems related with the entropic-force cosmology, see the recent [43| . 



20 



References 

[1] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1; 

V. Sahni, A. Starobinsky, Int. J. of Mod. Phys. A9 (2000) 373 | arX iv:astro-ph/9904398l ; 
T. Padmanabhan, Phys. Rep. 380 (2003) 235 [arXiv:hep-th/0212290] . 

[2] For a review, see e.g. P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75 (2003) 559 
|arXiv:astro- ph/0207347l ; 

E. J. Copeland, M. Sami, S. Tsujikawa, Int. J. of Mod. Phys. D15 (2006) 1753 
|arXiv:hep-th/0603057]; 

V. Sahni, A. Starobinsky, Int. J. Mod. Phys. D15 (2006) 2105 [arXiv:astro-ph/0610026 1 ; 
M. Li, X.-D. Li, S. Wang, Y. Wang, Commun. Theor. Phys. 56 (2011) 525 |arXiv:1103.5870] : 
D. Weinberg, M. Mortonson, D. Eisenstein, C. Hirata, A. Riess, E. Rozo [ ar Xiv: 1201. 2434 j . 

[3] B. Boisseau, G. Esposito-Farese, D. Polarski, A. A. Starobinsky, Phys. Rev. Lett. 85 (2000) 
2236 | arXiv:gr-qc/0001066j . 

[4] S. Nojiri and S.D. Odintsov, eConf C0602061 (2006) 06 [Int. J. Geom. Meth. Mod. Phys. 4 
(2007) 115] |arXiv:hep-th/0601213 |; 

T.P. Sotiriou, V. Faraoni, Rev. Mod. Phys. 82 (2010) 451, arXiv:0805. 1726 [gr-qc]; 

R. Woodard, Lect. Notes Phys. 720 (2007) 403 |arXiv:astr o-ph/0601672j ; 

S. Capozzieho, M. De Laurentis, Phys.Rept. 509 (2011) 167 [ arXiv: 1108.6266] . 

[5] L. Amendola, R. Gannouji, D. Polarski, S. Tsujikawa, Phys. Rev. D75 (2007) 083504 
|arXiv:gr-qc/0612180] ; 

W. Hu, I. Sawicki, Phys.Rev. D76 (2007) 064004 |arXiv:0705.lT58] : 

A. A. Starobinsky, JETP Lett. 86 (2007) 157 tarXiv:0706j04l]; 

S. A. Appleby, R. A. Battye, Phys.Lett. B654 (2007) 7 [a rXiv:0705.3T99] . 

[6] F. Bauer, J. Sola, H. Stefancic, JCAP 1012 (2010) 029 [a rXiv: 1006. 3944] : 
Phys. Lett. B688 (2010) 269 [ax Xiv:0912.0677] : 
Mod. Phys. Lett. A26 (2011) 2559 [arXiv: 1105. 1030] : 
Phys. Lett. B678 (2009) 427 [arXiv:0902.2215]; 

F. Bauer, Class. Quant. Grav. 27 (2010) 055001 [arXiv: 0909^237); 
Class. Quant. Grav. 28 (2011) 225019 [ar Xiv:1108.0875j : 

J. Sola, Fortsch. Phys. 59 (2011) 1108; 

S. Basilakos, F. Bauer, and J. Sola, JCAP 01 (2012) 050 [arXiv: 1109.4739] . 

[7] J. Sola, J. Phys. Conf. Ser. 283 (2011) 012033 [arXiv:1102.18T5] . 

[8] J. Sola, J. of Phys. A41 (2008) 164066 |arXiv:0710.4"T5T] . 

[9] I. L. Shapiro, and J. Sola, Phys. Lett. B682 (2009) 105 [arXiv:0910.4925j : 
confer also the extended version ,arXiv:0808.0315^ and references therein. 



21 



[10] O. Bertolami, Nuovo Cimento 93B (1986) 36; 

M. Ozer M. and O. Taha, Nucl. Phys., B287 (1987) 776; 
O. K. Freese K., et al., Nucl. Phys., 287 (1987) 797; 

J. C. Carvalho, J. A. S. Lima and I. Waga, Phys. Rev. D46, 2404, (1992); 

See e.g. the reviews: J. M. Overduin and F. I. Cooperstock, Phys. Rev. D58 (1998) 043506; 

and R.G. Vishwakarma, Class. Quant. Grav. 18 (2001) 1159, and references therein. 

[11] S. Basilakos, M. Plionis and J. Sola, Phys. Rev. D80 (2009) 3511 |arXiv:0907.4555] . 

[12] J. Grande, J. Sola, S. Basilakos, and M. Plionis, JCAP 08 (2011) 007 [arXiv: 1103.4632] . 

[13] J. Grande, J. Sola and H. Stefancic, JCAP 08 (2006) Oil |arXiv:gr-q c/0604057]; 
Phys. Lett. B645 (2007) 236 [arXiv:gr-qc/0609083j; 
J. Phys. A 40 (2007) 6787 [arXiv:gr-qc/0701090j . 

[14] E.P. Verlinde, JHEP 04 (2011) 029 [arXiv: 100 1.0785] . 



[15] T. Jacobson, Phys. Rev. Lett. 75 (1995) 1260 [arXiv:gr-qc/9504004] . 

[16] For a review, see e.g. T. Padmanabhan, Rept. Prog. Phys. 73 (2010) 046901 |arXiv:0911.5004] : 
Phys. Rep. 406 (2005) 49 ^arXiv : gr-qc /0311036j , and references therein. 

[17] D. A. Easson, P. H. Frampton, G.F. Smoot, Phys. Lett. B696 (2011) 273 [arXiv: 1002. 4278] . 

[18] D. A. Easson, P. H. Frampton, G.F. Smoot, Entropic Inflation [ arXiv: 1003. 1528] . 

[19] R. Casadio, A. Gruppuso, Phys. Rev. D84 (2011) 023503 iarXiv:1005.0790] . 

[20] T. S. Koivisto, D. F. Mota, M. Zumalacarregui, JCAP 02 (2011) 027 [arXiv:1011.2226] . 

[21] Y-F. Cai, J. Liu, H. Li, Phys. Lett. B690 (2010) 213 [arXiv: 1003.4526] : 
H. Wei Phys. Lett. B692 (2010) 167 [arXiv:1005JJ45]; 
Y.S. Myung, Astrophys. Space Sci. 335 (2011) 553 [arXi v: 1005.2240] : 
Y. Fu Cai, E. N. Saridakis, Phys. Lett. B697 (2011) 280 [arXiv: 101 1.1245] : 
M. Visser, JHEP 10 (2011) 140 [arXiv:1108.5240] . 

[22] I. L. Shapiro, and J. Sola, JHEP 02 (2002) 006 [arXiv:he p-th/0012227] ; 
Phys. Lett. B475 (2000) 236 [arXiv:hep-ph/9910462]; 
Nucl. Phys. Proc. Suppl. 127 (2004) 71 [arXiv:hep-ph/0305279] ; 
JHEP proc. AHEP2003/013 [arXiv:astro-ph/0401015l; 



J. Sola, H. Stefancic, Phys. Lett. B 624 (2005) 147 |arXiv:astro-ph/0505133j ; 
Mod. Phys. Lett. A21 (2006) 479 [arXiv:astro-ph/0507110 |. 

[23] S.W. Hawking, G. T. Horowitz, Class. Quant. Grav. 13 (1996) 1487 |arXiv:gr-qc/9501014] . 

[24] T. Padmanabhan, Class. Quant. Grav. 22 (2005) L107 [ arXiv:hep-th/0406060] . 



22 



[25] M. Maggiore, Phys. Rev. D83 (2011) 063514 |arXiv:1004.rf82] : 
N. Bilic, Phys. Rev. D83 (2011) 105003 [arXiv: 1104. 1349] : 
L. Hollenstein, M. Jaccard, M. Maggiore, E. Mitsou, jarXiv:1111^5575] : 
Phys. Lett. B704 (2011) 102 ^arXiv:1104. 3797j: 

N. Bihc, S. Domazet, B. Guberina, Phys. Lett. B707 (2012) 221 [ arXiv: 11 10.2393] . 
[26] F. E. M. Costa, J. A. S. Lima, F. A. Ohveira (arXiv:1204.1864j . 

[27] H. Fritzsch, and J. Sola, Matter non-conservation in the Universe and dynamical dark energy 
[arXiv:1202.5097j. 

[28] K. Jedamzik, M. Pospelov, New J. Phys. 11 (2009) 105028 [arXiv:0906.2087] . 
[29] R. Amanullah et al., Astrophys. J., 716 (2010) 712 |arXiv:1004.mT] ]. 
[30] S. Nesseris and L. Perivolaropoulos Phys. Rev. D. 72 (2005) 123519 
[31] H. Wei and Rong-Gen Cai Phys. Lett. B. 663 (2008) 1 

[32] D. J. Eisenstein et al. (SDSS Collab.), Astrophys. J., 633 (2005) 560 | arXiv:astro-ph/0501171] ; 
N. Padmanabhan, et al. (SDSS Collab.), Mon. Not. Roy. Astron. Soc, 378 (2007) 852 
[a?Xivrastro-ph/0605302| . 

[33] W. J. Percival et al. (SDSS Collaboration), Mon. Not. Roy. Astron. Soc, 401 (2010) 2148 
[arXiv:0907.1660J; 

E. A. Kazin, Astrophys. J., 710, 1444 (2010) t arXiv:09 08.2598j. 

[34] C. Blake et al., Mon. Not. Roy. Astron. Soc., 418 (2011) 1707 [ar Xiv:1108.2635] . 

[35] D. J. Eisenstein, and W. Hu, Astrophys. J. 496 (1998) 605 (arXiv:astro-ph/9709112| . 

[36] J. R. Bond, G. Efstathiou and M. Tegmark, Mon. Not. Roy. Astron. Soc. 291 (1997) L33 
[a rXiv:astro-ph/9702100| . 

[37] S. Nesseris and L. Perivolaropoulos, JCAP 0701 (2007) 018 [ arXiv:astro-ph/0610092] . 

[38] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Supl. 180 (2009) 330 
[arXi v:0803.0547j : 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Supl. 192 (2011) 18 
[arXi v: 1001. 4538) . 

[39] O. Elgaroy, T. Multamaki, Astron. Astrophys. 471 (2007) 65 [arXiv:astro-ph/0702343] . 

[40] T. M. Davis et al. Astrophys. J. 666 (2007) 716 

[41] M. Plionis et al. Mon. Not. Roy. Astron. Soc. 416 (2011) 2981 

[42] M. Hicken et al., Astroplys. J., 700 (2009) 1097 [arXiv:0901.48041 . 

[43] J. J. Roveto, G. Munoz [a rXiv:1201.2475] . 



23 



[44] J. Grande, A. Pelinson, and J. Sola, Phys. Rev. D79 (2009) 043006 [arXiv:0809.3462] : 



24 



